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Abstract. We present a conjecture on the exact location of the multicritical point 
in the phase diagram of spin glass models in finite dimensions. By generalizing our 
previous work, we combine duality and gauge symmetry for replicated random systems 
to derive formulas which make it possible to understand all the relevant available 
numerical results in a unified way. The method applies to non-self-dual lattices as well 
as to self dual cases, in the former case of which we derive a relation for a pair of values 
of multicritical points for mutually dual lattices. The examples include the ± J and 
Gaussian Ising spin glasses on the square, hexagonal and triangular lattices, the Potts 
and Z q models with chiral randomness on these lattices, and the three-dimensional ± J 
Ising spin glass and the random plaquette gauge model. 



1. Introduction 

Properties of finite-dimensional spin glasses are still under debate although the problem 
is essentially settled for the mean-field model pQ. Outstanding problems for finite- 
dimensional spin glasses include the existence or absence of spin glass phase and whether 
or not the mean-field picture of the spin glass phase applies. Another interesting, but 
less extensively studied, issue is the structure of the phase diagram, in particular where 
precisely the multicritical point is located and what the values are for the critical 
exponents characterizing the system behaviour at and away the multicritical point. 
The present paper discusses this problem of the location of the multicritical point for 
finite-dimensional spin glass models by analytical methods. 

A number of numerical investigations on this problem exist for various lattices. 
However, it has been quite difficult to derive analytical results for regular finite- 
dimensional lattices until a few years ago when we succeeded in devising a method 
to predict the exact locations of the multicritical points for the square lattice Ising and 
Potts models and four-dimensional random plaquette gauge model using duality, gauge 
symmetry and the replica method El E] I n the present paper we generalize this 
theoretical framework so that it is applicable to a pair of mutually-dual lattices, for 
which our theory relates the pair of values of multicritical points of the two lattices. 

The logic of our theory includes a step which is yet to be justified rigorously, and 
hence the status of our result is a conjecture at this moment. Nevertheless, our theory 
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enables us to understand all the relevant available numerical data for the multicritical 
points derived independently by a number of authors. Also our theory satisfies necessary 
conditions which the exact solution should satisfy. 

This paper is organized as follows. In the next section, we recall the basic 
formulation which was developed in our previous studies in order to fix the notation 
and set the stage for further developments in the following sections. In section 3, 
a generalization of the theory to non-self-dual Z 2 (Ising) models with randomness is 
discussed. This argument is followed by section 4, in which we further generalize 
the arguments to Z q systems with chiral randomness. The final section is devoted 
to conclusion and discussions. 



2. Self-dual Z 2 models 

In this section we briefly review the duality arguments in [3J H] applied to the two- 
dimensional (2d) ± J random bond Ising model on the square lattice and the 4d random 
Z 2 lattice gauge model, by which we fix the notation and set a stage for generalizations 
in the following sections. After recalling the duality arguments for non-random systems, 
we apply the idea to random systems. 



2.1. Duality of non-random models 

We first elucidate the duality of a non-random Z 2 (Ising) system [Hj. Let us prepare a 
c?-dimensional lattice and assign Z 2 spins on r — 1 dimensional elements x on the lattice, 
which we denote by Sx- We consider a model on the lattice whose Hamiltonian is given 
by 

C X£dC 

where C is the r-dimensional element on the lattice and dC is its boundary of dimension 
r — 1 5|. Let u±i denote the Boltzmann factor for an element C, 

u ±1 (K)=e ±K , for J] S x = ±1, (2) 

xedc 

where K = [3J. For the case of the usual Ising model (r = 1), U\{K) is the bond (edge) 
Boltzmann factor for parallel spins at both ends and U-i(K) is for anti-parallel spins. 
Then the partition function Z is a function of u±\\ Z = Z{u\(K),U-i(K)}. 

The dual model is defined on the dual lattice (the definition of which is given in j^j). 
The dual Boltzmann factor for the dual element C* is defined by the two-component 
Fourier transformation of u±\ 

, ,,^_ u l {K)±U- 1 (K) e K ± e- K 
u ±1 (K) = ^ = (3) 

For the present Hamiltonian (|T|). the dual Hamiltonian is given by 

h*=-jyi n s *> w 

c* xedc* 
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which is of the same form as equation (|T|). C* is the dual element of C and has dimension 
d — r. Next we derive the dual expression of the partition function for general dual pairs. 
The partition function as a function of u has the following property, 

Z olig { Ul (K),u^(K)} = 2 a Z dnal {ul(K),u*_ 1 (K)}, (5) 

Here Z orig is the partition function of the original model and Z dual is for the dual model. 
The symbol a is a constant determined by the numbers of the elements on the lattice. 
(See appendix.) 

We give several examples of dual pairs below, which would be helpful to understand 
the definition of the dual transformation. The arrows indicate duality relations. 

• (d — 2, r — 1) 2d Ising model on the square lattice <-> 2d Ising model on the square 
lattice (self-dual) 

• (d — 2, r — 1) 2d Ising model on the triangular lattice 2d Ising model on the 
hexagonal lattice 

• (d = 4, r = 2) 4d Z 2 lattice gauge model on the hypercubic lattice «-> 4c? Z 2 lattice 
gauge model on the hypercubic lattice (self-dual) 

• (d — 3, r = 1) 3d Ising model on the cubic lattice <-» 3d Z 2 lattice gauge model on 
the cubic lattice 

So far the discussions have not been restricted to self-dual models. 

For self-dual cases, Z OTig and Z dual are the same function and the prefactor 2 a 
on the right hand side in (JHJ) becomes a trivial constant which is negligible in the 
thermodynamic limit (and is omitted in the following). Hence (0) is simplified to 

Z{ Ul (K),u^(K)} = Z{vZ(K),u*-i(K)h (6) 

where the symbol Z or ± g / dua i is simplified to Z. From this expression it is seen that Z is 
invariant under the exchange U\(K) <-> u\(K) and u_i(K) <->• u^^K), which means self- 
duality of the partition function. The critical point of a self-dual model is obtained, if it 
is unique, by the fixed-point condition for these Boltzmann factors, u±i(K c ) = w^_ 1 (if c ), 
which yields K c = ~ ln(v2+ 1). It is clear that this transition point is shared by the 2d 
Ising model on the square lattice, the 4d Z 2 gauge model on the hypercubic lattice and 
their higher dimensional generalizations [S]. 

2.2. Duality of random models and conjecture for the critical point 

Let us introduce randomness. To investigate critical points of random systems with the 
aid of the duality formalism, we will utilize the technique of reference |2] with some 
modifications. 

The random model treated here is a system with Z 2 variables and bimodal 
randomness. The Hamiltonian is written as 

h=-jj2 t c n s *> ( ? ) 

C XedC 
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where Tc is a quenched random variable dependent on each element C. r c takes the 
value 1 with probability p and —1 with 1 — p. To treat random systems, we employ 
the standard replica method. Let us consider the n-replicated system and define the 
averaged Boltzmann factor Xk for an element C, which corresponds to the configuration 
ricEeac &x = 1 in n — A; replicas and — 1 in k replicas. The explicit form of Xk is 

x k (p, K) = pe {n ~ 2k)K + (1 - p)e^ n - 2k)K . (8) 

The n-replicated partition function is, after average over randomness, a function of these 
Boltzmann factors, 

[Z n ] w = Z n {x (p, K),x\(p, K),..., x n (p, K)}, (9) 

where [ ] av means random average. 

We also define the dual Boltzmann factor x* k (p, K) on the dual lattice. The explicit 
forms are obtained by the two-component Fourier transformation with the result 

x* 2k (p, K) = 2- n '\e K + e- K ) n - 2k (e K - e~ K f\ 

x* k+1 (p, K) = 2- n ' 2 (2p - l)(e K + e- K ) n - 2k -\e K - e - K ) 2k+1 , (10) 

where k is a non-negative integer in the range < 2k < 2k + 1 < n. The partition 
function satisfies a generalization of 

where a is an appropriate constant. Now we restrict our attention to the case where the 
system is self-dual when we remove randomness (e.g. 2c? ±J random bond Ising model 
on the square lattice). Using (fTT]) we can express the self duality of the n-replicated 
partition function for such a system as 

Zn{_%0i • • • i -^n} = ^"n\_%Qi • • • > ^nf > (12) 

where an overall constant is neglected. Thus self duality is recognized by the fact 
that Z n is invariant when the exchanges Xk{p,K) <-> x* k (p,K) for all k are performed 
simultaneously. 

It is in general impossible to identify the transition point from the fixed-point 
condition of the duality relation ([12)1 . unlike the non-random case because the fixed- 
point conditions of all the variables not satisfied 
simultaneously. The authors of j3 Hj nevertheless developed an argument leading to a 
conjecture that the fixed-point condition of the leading Boltzmann factor 

x (p c , K c ) = x* (pc, K c ) (13) 

can well be the most plausible candidate to give the exact transition point of the random 
system at least on the Nishimori line (NL) [Zj, e~ 2K = (1 —p)/p, where enhanced 
symmetry simplifies the system properties significantly. 

This prediction has been confirmed to be correct in the cases of n = 1,2 and oo jHj- 
It has also been shown that numerical results for the quenched limit n — > |Hj agree 
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very well with the conjectured value of p c = 0.889972... which is the solution to the 
formula obtained in the n — > limit of Xq = x* on the NL: 

-plog 2 p- (l-p)log 2 (l-p) = -. (14) 

This conjecture also leads us to an interesting result that the multicritical points of the 
2d ±J random bond Ising model and the Ad random Z 2 lattice gauge model (random 
plaquette gauge model) are located at the same point on the p-K plane This 
observation has also been confirmed numerically by a recent study jH]- 

Determination of the multicritical point is also quite important from the standpoint 
of the quantum information theory, quantum memory in particular. To be specific, the 
value 1 — p c at the multicritical point of the 2c? random bond Ising model is equivalent to 
the accuracy threshold of the 2d toric code with perfect measurement [TU1 lllj. which is 
estimated to be 0.110028 . . . from the above discussion. In addition, the value 1 — p c at 
the multicritical point of the 4d random gauge model also gives the accuracy threshold of 
the 4d toric code (or 3d code with imperfect measurement) [I], which is also determined 
as 0.110028.... 

3. Z 2 models 

With the knowledge of the previous section in mind, we proceed to discussions on the 
non-self-dual cases with and without randomness. 

3.1. Duality for non-self-dual cases 

In this subsection, we develop an argument for duality of generic non-self-dual Z 2 models. 
It is clear that a simple duality relation like (jHJ) is not enough to determine the transition 
point. Nevertheless we show that, by generalizing the arguments in the previous section, 
one can still derive a relation between the transition points of a model and its dual. 

First let us discuss non-random systems. Consider the product of partition 
functions of the original and dual models with inverse temperatures K\ and K 2 , 
respectively. From one finds 

= Z olig {ul(K 2 ),uU(K 2 )}Z A uA<(Ki),u*_ 1 (K 1 )}, (15) 

which indicates that the product is invariant under the simultaneous exchange 
u±i(Ki) <-> u* ±1 (K 2 ) and u±i{K 2 ) <-> u* ±l {K\). Hence, if there is a unique transition 
point Ki c (resp. K 2c ) in the original (resp. dual) model, the relation between two 
critical points is expected to be given by 

u ±1 {K lc )u ±1 {K 2c ) =u* ±1 {K lc )u* ±1 {K 2c ), (16) 

which is invariant under the transformation above. One can verify that this is equivalent 
to e~ 2K2c = tanh K\ Cl which gives the correct relation between the two transition points. 
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Now we move on to random systems. In this case, we can express the duality 
between two random models (fTTj) as follows, 

2n,orig{£o(Pl, Kl), ■ ■ ■ , Xn(Pl, K^} Z n:du£k \{x (p 2 , K 2 ),..., X n (p 2 , K 2 )} 

= Z niOrig {x* (p 2 , K 2 ),..., x* n (p 2 , K 2 )}Z n:dual {x* (p u K x ),..., x* n (p u Kx)}, (17) 

where p\ , p 2 and K\ , K 2 denote the probability of positive interaction and the inverse 
temperature for the original/dual models, respectively. Thus the product of the partition 
functions is invariant under the simultaneous exchange Xk(pi,Ki) x* k (p 2 ,K 2 ) and 
Xk(P2, K 2 ) «-> a£(pi, K x ) for all k. 

The argument developed so far naturally suggests that the relation between the 
critical points of the original and dual systems is given by the fixed-point condition of 
the leading Boltzmann factors at least on the NL. Explicitly, this condition reads 

x (pic, K lc )x (p 2c , K 2c ) = x* (p lc , K lc )x* Q (p 2c , K 2c ), (18) 

in conjunction with the NL condition 

e- 2 ^ = e" 2 ^ = i^. (19) 

Pi P2 

Equation ()18|) with ()19|) is written in terms of p\ c and p 2c as 

( Plc n+1 + (1 - j9 lc )™ +1 ) (p 2c n+1 + (1 - p 2c ) n+l ) = 2- n . (20) 

We expect that the two multicritical points are related by this equation. If the quenched 
(n — > 0) limit is taken, this yields the relation, 

H(p lc ) + H(p 2c ) = 1, (21) 

where 

H(p) = -plog 2 p- (1 -p)log 2 (l -p). (22) 

Therefore our main statement for the non-self-dual Z 2 model is the following: the two 
critical points on the NLs of the mutually-dual systems with quenched randomness are 
expected to be related by equation (|2"T]) . 

There are many reasons to believe that our conjecture expressed in f!20|) and (|2*T|) 
is exact. Some of them are given in the rest of this paper. For simplicity the discussion 
in the following is restricted to the 2d random bond Ising models on the mutually-dual 
lattices such as the hexagonal and triangular lattices, though the conjecture applies 
quite generally to arbitrary systems described by the Hamiltonian (JZJ). 

Before closing the present subsection, we explain an explicit representation of the 
dual random models [21 E], which is necessary in the following discussions. The ratios 
of dual Boltzmann factors (fTU|) to Xq are 

x *2k/ x o = tanh 2fc if, 

4fe+i/4 = ( 2 P - 1) tanh 2fc+1 K. (23) 

These ratios of Boltzmann factors are realized by a system with the following explicit 
Boltzmann factors written in terms of Ising spin variables 

Aexp (K(SW + S (2) + . . . + S (n) ) + K P S W SW . . . , (24) 
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where is the product of Ising spin variables in the kth replica. For example, in the 
case of the usual nearest neighbour interactions, stands for s!- Sj . K and K p are 
defined by 

e~ 2k = tanh K, e~ 2Rp = 2p - 1 (= tanh K p ) . (25) 

From the expression ()24j) . the dual of the ±J random bond Ising model can be 
interpreted as the model with the non-random Ising interaction in each replica and 
the interaction between replicas. If the condition of the NL, K = K p , is imposed, this 
turns to 

Aexp (k(SW + S (2) + . . . + S {n) + S {1) S {2) . . . S (n) )) . (26) 
3.2. Verification for n = 1 

Let us first show that the relation (fT%j) gives the exact answer even without the NL 
condition when n — 1. Equation (fTHj) is, for n — 1, 

( Pc e K ^ + (1 -pi c )e- K ^)(p 2c e K *' + (1 -J9 2c )e~^) = ~(e^ + e~^)(e x - + e~^).(27) 



^From the expression (j24j) . it is found that the ±J random bond Ising model for n = 1, 
averaged over randomness, is regarded as a non-random Ising model with coupling 
K + K p on the dual lattice. For example, the random bond Ising model on the 2d 
triangular lattice with parameters pi and Ki is equivalent to the non-random Ising 
model on the hexagonal lattice, whose inverse temperature K\ is given by 

K x = k x + K vU (28) 

where K v \ is defined by the second expression of (f23|) with p replaced by pi . Conversely, 
the random bond Ising model on the 2d hexagonal lattice with parameters p 2 and K 2 
corresponds to the non-random Ising model on the triangular lattice, whose inverse 
temperature K 2 is 

K 2 = K 2 + K p2 . (29) 

As is well-known, the non-random Ising models on the 2d triangular and the hexagonal 
lattices are mutually- dual, and two critical inverse temperatures satisfy [T2] 

e~ 2klc = tanhK 2c . (30) 

Using equation (}25|) we can confirm the equivalence between equations (|27]l and (J30|) . 
Note that the condition of the NL is not used, and this equivalence holds everywhere 
on the phase boundary. 

3. 3. Verification for n = 2 

Next we discuss the two-replica case. Equation ()18|) for n = 2 is 

(pi c e 2 ^ + (l-pi c )e- 2X -)(P2 C e 2 ^ + (l-p 2c )e- 2X2c ) = i( e ^+ e -^) 2 ( e ^+ e -^) 2 .(31) 
In the rest of this subsection, we restrict our attention onto the NL. 



Exact location of the multicritical point: A conjecture 



8 



Due to the dual representation (J26)) for n = 2, we know that the random bond 
Ising model is equivalent to the non-random four-state Potts model on the dual lattice; 
four states are constructed by the direct product of Ising factors in two replicas 
S« = ±1,S® = ±1, and the dual Boltzmann factor is Ae zk for = = 1 
and Ae~ K otherwise. Therefore the system corresponds to the four-state Potts model 
on the dual lattice with coupling 2K, namely, 

(3H = -2K (25 5(1 ) il V) )1 - 1) - K. (32) 

For example, the random bond Ising model on the triangular lattice is equivalent to 
the non-random four-state Potts model on the hexagonal lattice and vice versa. Here 
we denote the parameters of the random bond Ising model on the triangular lattice by 
Pi,K\ and that on the hexagonal lattice by p2,K2. Then the inverse temperatures of 
the corresponding four-state Potts models Ki,K 2 are given by 

K x = K 2 = 2K 2 , (33) 

as mentioned above. The dual coupling K\ (or K 2 ) is defined by equation (f2"5j) . 

The inverse critical temperatures of the four-state Potts models on mutually-dual 
lattices satisfy [12J 

e -2^ 2c = ± (34) 

e K ^ + 3e~ K ^ 

It is straightforward to show that equation (}3*Tj) is equivalent to equation (J53j) under the 
NL condition (fTTl. 



3.4- The limit n — ► oo 

Next is an argument for the n — > oo limit. The partition function for arbitrary n is 
expressed as 

Zn = [Z n U = [e- n0Nf{K) U (35) 

where N is the number of sites and f(K) is the free energy for the inverse temperature 
K. If the replica number n is taken to be quite large, the averaged partition function is 
dominated by the configuration of the smallest value of f(K), which is realized by the 
perfect ferromagnetic bond configuration and its gauge equivalents. Using this argument 
Z n can be approximated as 

2 n «e- n »W, (36) 

where fo{K) is the free energy of the non-random ferromagnetic system. Therefore the 
random bond Ising system for n — > oo is interpreted as the non-random ferromagnetic 
Ising system (and its critical point is denoted by Ki c ). 

We also consider the random bond Ising model on the dual lattice, and for n — > oo 
we can regard it as the non-random Ising model with the critical point K 2c . These two 
critical points are related by the well-known equation, 

e - 2K ^ = tardive. (37) 
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Figure 1. Relaxation of magnetization near the multicritical points of the 
random bond Ising model on the hexagonal (left) and the triangular (right) 
lattices from NER analysis. We prepared L x L(L = 10 3 ) spins and averaged 
the results over 200 samples. The system at criticality is expected to yield a 
straight line in this log-log plot. 



If we consider the n — > oo limit of equation (fTHj) combined with the NL condition ([19)1 . 
we obtain 



D K lc +K 2c 



+ e 



){e K ^ + e- K ^) 



(3f 



which is equivalent to equation (1371) . 



3.5. Numerical evidence 

We have checked the relation (|2*T|) in the quenched limit numerically. We executed 
Monte Carlo simulations for the 2d random bond Ising models on the hexagonal and 
the triangular lattices, a dual pair. To observe the criticality, we make use of the non- 
equilibrium relaxation (NER) method which yields the power-law behaviour of 
decreasing magnetization with Monte Carlo steps on criticality. For this method we 
prepare all-up spins as the initial state and let them relax in each Monte Carlo step. In 
order to identify the multicritical point, we choose the parameters of the system to be 
on the NL and vary the parameter p (and K is also varied accordingly). 

The results are shown in figure 1. From these results, the locations of the 
multicritical points are estimated at pi c = 0.930(5) (or 0.347 < H(px c ) < 0.384) for the 
hexagonal lattice and at p 2c = 0.835(5) (or 0.634 < H(p 2c ) < 0.658) for the triangular 
lattice. The two H's sum up to 0.981 < H(pi c ) + H(p 2c ) < 1.042, a consistent result 
with our conjecture in equation (|2"Tf . 



3.6. Random models in 3d 

As already noted in section 2, the random bond Ising model and the random Z 2 lattice 
gauge model on the 3d cubic lattice are mutually dual. Thus we expect that the 
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multicritical points of these two models will satisfy equation (f2T|). This conjecture 
can be confirmed by numerical simulations of these 3d systems. 

In fact the multicritical points of these two models have already been estimated. For 
the 3d random bond Ising model the multicritical point is estimated to be p\ c = 0.7673(3) 
[Hj which yields H(pi c ) 0.783. For the 3d random gauge model it is p2 C ~ 0.967 [THj . 
giving H(p 2c ) ~ 0.209. From these results the sum H{pi c ) + H{p 2c ) is about 0.992, a 
reasonable value in view of our expectation (}2"Tj) . 

From the viewpoint of the quantum information theory, the accuracy threshold of 
the 2d toric code with imperfect measurement can be determined by the value 1 — p c 
at the multicritical point of the 3d random gauge model ^13 ITT1 . Therefore it is an 
advantage of the present analysis that the duality gives an alternative way to determine 
the accuracy threshold when we utilize the result of the 3d random bond Ising model. 
From the result of the direct numerical analysis of the 3d random gauge model ^Hj and 
the duality argument combined with the numerical result of the 3d random bond Ising 
model [Hj, the accuracy threshold is estimated to be about 0.03 ~ 0.035. 

3. 7. 2d anisotropic random bond Ising model and self duality 

In this and the next subsections we discuss the 2d random bond Ising model with 
anisotropic disorder. This model is not a non-self-dual model, but the structure of the 
problem is similar to the mutually-dual case. We show that the duality formalism can 
be applied to this system and gives a conjecture for the critical points. It should also 
be noted that the critical point of this model gives the accuracy threshold of the Id 
quantum repetition code with imperfect measurement [TT] . 

First we review the duality of the anisotropic non-random Ising model on the 2c? 
square lattice. The Hamiltonian is 

H = —Jh SiSj — J v ^ SiSj, (39) 

where Jh and J v are uniform coupling constants. The symbol Ch denotes the set of 
horizontal bonds and C v is for vertical ones. The partition function is 

z = e n n ( 4 °) 

{Si=±l} (ij)EC h (ij)ec v 

where Kh V = f3Jh, v — Jh,v/kT and Sij = S^Sj. The symbol u is the Boltzmann factor 
for the bond between i and j, 

u 3 ..(K hjV ) = exp(K htV S i:j ). (41) 

Next we define the dual Boltzmann factor by the binary Fourier transformation, 

^_ u 1 (K h>v )±U- 1 (K htV ) e K ^ ± e~ K ^ 
u ±l {K Kv ) = -j= = -j= . (42) 

Using these Boltzmann factors, we can express the duality of partition function as 

Z{u 1 {K h ),u-. 1 (K h ),u 1 (Kv),u- X (K v )} = Zi^K^ut^K^iK^ut^Kh)}. (43) 
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In this expression the Boltzmann factors of the vertical and horizontal bonds are 
exchanged because a vertical bond is mapped to a horizontal bond on the dual lattice 
and vice versa. 

It is obvious that the partition function is invariant under the exchange u±\(Kh) <-» 
u±i(K v ) and u±i(K v ) <->• u± x {Kh), which is similar to the case of mutually-dual non- 
random systems. The critical points are determined by the equation, 

u ±1 {K h )u ±l {K v ) = u* ±1 {K v )u* ±1 {K h ), (44) 

which yields e~ 2Kv = tanhi^. 

Next we study the system with randomness. The Hamiltonian is 

H = -J v r^SiSj-Jn ]T r^S % S v (45) 

{ij)eC v {ij)&C h 

where T v,h are random variables which depend on each bond and obey the probability 
distribution, 

Ptft v ) = Ph,A^ v - 1) + (1 - Pn,v)5(4' v + 1), (46) 
for horizontal (h) or vertical(f) bond. The averaged partition function is a function of 
averaged Boltzmann factors x k (ph, K h ) and x k (p v , K v ), 

[2"% v = Z n {x (p h ,K h ), . . . ,x n (p h ,K h ),x (p v ,K v ), . . . ,x n (p v ,K v )}. (47) 
We also define the dual averaged Boltzmann factors in the same way as in section 2, 

x* 2k {p h , v , K h , v ) = 2- n ' 2 (e K ^ + e - K ^) n - 2 \e K ^ - e~ K ^) 2 \ 

X* 2k+1 ( Ph , v , K h>v ) = 2~ n / 2 {2 PKv - l)(e K ^ + e -K h , v) n-2k-X {e K h , v _ e -K h , v) 2k + ^ (4g) 

Using Xk and x* k , we can express the duality of the n-replicated partition function, 

Z n {x (ph, K h ), • • • , x n (p h , K h ), x (p v , K v ), x n (p v , K v )} 

= Z n {x* (p v , K v ),..., x* n (p v , K v ),x* (p h , K h ),..., x* n (p h , K h )}, (49) 

which is invariant under the simultaneous exchange x k {ph,K h ) <-> x* k (p v ,K v ) and 
x k (p v , K v ) x* k {ph, Kh) for all k. The Boltzmann factors of the vertical and horizontal 
bonds should be exchanged as in the non-random case. 

From the argument above, we make a conjecture for the critical points on the NLs, 

x {p h , K h )x (pv, K v ) = x* (p h , K h )x* (p v , K v ), (50) 

from the analogy with the non-random case f)44j) or the non-self-dual random case (fTSj) . 
In this system, the NLs are defined for the horizontal and vertical bonds, respectively, 
by 

e -2K* = l^Ph^ e -2K* = l^Pv (51) 
Ph Pv 

We should note that the original system has four parameters ph, K h ,p v , K v and the two 
conditions of the NLs reduce the number of independent variables to two. Consequently, 
the conjecture f)50|) combined with the conditions of the NLs (JSTj) does not fix the 
multicritical point on the phase diagram but determines the location of the "critical 
line". This fact is favourable for numerical verification of our conjecture because one 
parameter can be chosen freely even on criticality and on the NLs. 
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3. 8. Verification for anisotropic system 

We verify the validity of the conjecture in equation (}50|) . For this purpose the following 
property is useful; the form of equation ()5U)) is completely the same as equation (fTH|) 
when the parameters are replaced as 



Using this correspondence we can check the conjecture ()50|) by the same argument as 
in mutually-dual systems with n = 1, 2 and n — > oo. 

• n = 1: The anisotropic random bond Ising model is equivalent to the anisotropic 
non-random Ising model, and equation (|30)l turns to the critical condition for the 
anisotropic non-random model when the replacement is used. 

• n = 2: The system is equivalent to the anisotropic non-random four-state Potts 
model. Using the replacement (JH2J), equation (|3*4*)l becomes the critical condition 
for the anisotropic Potts model. 

• n — > oo: The system is equivalent to the anisotropic non-random Ising model and 
the critical point satisfies equation (|3*7j) using the correspondence (152)) . 

If we consider the quenched limit n — > 0, we obtain the relation, 



by the same argument as in mutually-dual systems (j21)l . Here we take Ph,Pv as the two 
independent variables and eliminate K^ v using the conditions of the NLs. We expect 
that ph and p v will satisfy equation (f53|) on criticality. 

We have checked equation (|53*)) numerically. The result is shown in figure 2. 
The locations of the critical points expected from equation (153)1 are in reasonable 
agreement with numerical results. Using the NER with ordered and disordered initial 
configurations, we also checked that the result for p v very close to 1 (or p h ~ 1) does 
not show discrepancy from the conjecture, which is not shown in the figure because the 
critical point could not be determined with sufficient precision due to rigid spin domains 
along one axis. Thus we conclude that our conjecture is consistent with the numerical 
result for the anisotropic system as well. 

We are allowed to interpret (|53*)l as the accuracy threshold of the Id quantum 
repetition code, given by the value of 1 — ph under the imperfect measurement rate 
1 -p v . 

4. 2d Z q model 

In this section we consider a generalization of the formalism to a multi-valued spin 
model, the random Z q model. The Z q model has q states per site. In the non-random 
case, this system in two dimensions has both aspects of the ferromagnetic Ising model 
(i.e. with ferromagnetic order) and the XY model (i.e. with the Kosterlitz-Thouless 
(KT) order) for q not too small ^OlE)- The authors of reference [Tj\ showed that this 



{p 1 ,K 1 ,p 2 ,K 2 } -> {p h ,K h ,p v ,K v }. 



(52) 




(53) 
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Figure 2. Results by the NER for the anisotropic random bond Ising model. We 
simulated systems with up to 500 x 500 spins under all-up initial configuration. 
The results are averaged over 200 samples. For executing numerical calculation, 
we fixed p v and scanned p^ in order to search the critical point. The critical 
points obtained are shown with error bars. The curve H(ph) + H{p v ) = 1 lies 
within the error bars for all p v . 




model has three phases for q > 5, disordered, KT-like, and ferromagnetic, if there exists 
a KT phase in the continuum limit q — > oo. Duality has been used to relate the two 
transition points. 

The objective of the present section is to generalize the duality arguments in the 
previous section and apply them to the randomized version of the Z q model. As far as 
the authors know, this type of models have not been studied analytically in detail in 
the literature. 

4-1. Duality and analysis of the transition points 

We consider the random chiral Z q model, for which the partition function is written in 
the form, 

z = J2Y[e VK ^- k i- h >\ (54) 

where Vx(k) is the interaction which satisfies cyclic condition, Vxik + q) = V^k). Here 
i,j are the site indexes, k is the state variable which takes q values, 0, 1, • • • , q — 1, with 
q G Z + = {1, 2, • ■ ■} and I is the random variable on each bond which also takes q values, 
0, 1, • • • , q — 1. The probability that a random variable / on a bond takes the value I is 
denoted by pi (J2iPi = !)■ A collection of them is denoted by {p}(= {po,Pi ■ ■ -Pq-i}) 
in the sequel. The summation Yl{k} ^ s taken over state variables on all sites and the 
product Yliij) i s taken over all bonds. Only nearest neighbour interactions are assumed 
to exist on the square lattice in the present section. 
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Our discussions are mostly for the Z q version of the Villain model ^B|, for which 
the Boltzmann factor is given by 



e V K (k) 



e -f (7 fc " 2 ™) 2 ( = _L= \^ e -l 2 /(2K)+2nilk/q \ _ ^ 
m =-oo \ l& J 

We generalize the duality arguments to the random Z q model. Let us average the 
n-replicated partition function. The resulting function Z n can be written in terms of 
the q n variables (local Boltzmann factors for neighbouring spin pairs), 

Xk u -,k n ({p},K) = J2P^ VK(l+kl)+ '" +VKil+kn ^ ( 56 ) 
i 

where ki = 0, 1, — 1 (i = 1, 2, • • ■ , n). The variable Xfci,-,fcn({p}> K) is the 
generalization of the Boltzmann factor Xk(p,K) which appeared in section 2. If we 
set q = 2 and Vrt(0) = K, V^(l) = — K, we obtain the relation between the present and 
previous Boltzmann factors, 

n 

Xki,-,k n {{p}, K) = x k (p, K), where k = ^k m {k m = 0, 1). (57) 

m=l 

We can define the dual model, for which the partition function can be written in 
terms of the dual variables. They are defined by g-value Fourier transformations jH] 

Xk,., kn ({p}, K) = g"" /2 J> E ^ fcimi e^ (mi+i) • • ^e v ^ +l \ (58) 

l mi m n 

where to = e 2m ^ q . A remark is in order. The variables X&i,-,fcn are already defined at 
the end of section 2 of reference jS] in a different form from the ones given above. The 
expressions (74)-(79) in [3| are not in general valid and should be replaced by (JSTj) and 
(PI). 

Now we define 

Xo({p},K) = X o,., ({p},K) = J>e^« 



q-l /g-1 \ n 

X* ({P}, K) = X*a,..,o({p}, K) = <T" /2 $> £ e v *™ , (59) 

1=0 \k=0 J 

and apply the same arguments as in the previous section to identification of transition 
points of the Z q model. If uniqueness of the critical point is assumed, generalization of 
(fT3|) gives 

Xo({p},K)= X *o({p},K). (60) 

Remember that the non-random model does not have a unique transition point for 
q > 5 [17]. Rather, two transition points are related to each other by the duality 
relation. Therefore let us consider the case where there exist two transition points for 
the random case. If there are two critical points K\, K 2 for a given {p}, they may satisfy 

X o({p}, Kjxodp}, K 2 ) = x *o({p}, Ki)x*o({p}, K 2 ). (61) 
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Let us take the quenched (n — > 0) limit of these relations. When uniqueness is 
assumed (K\ = K 2 ), the condition reads 

£=i / e v *« \ 1 

"S P,1 ° g 'lEg^j = 2- (62) 

If there are two critical points, this is replaced by 

£i , / e^W eW) \ , . 

l^PitoE, ^ q -i e y Ki(k+l) ■ vg -i eWfc+ o = ( 63 ) 

If we consider the non-random case, where p — 1, = 0(Z 7^ 0), %o an d Xo re& d 
XO ({1,0...},^)= ^e-fCW, (64 ) 

X&({1,0 . . .}, K) = g ^ e-f <t"™>* = 1 £ e -#. (65) 

If there are two critical points, equation (|61|) gives the relation between them, 

K X K 2 = (66) 

This agrees with the correct relation However, it is in general not expected that the 
condition (JBTj) determines the whole shape of the phase boundary in the phase diagram 
of the random Z q model. Restriction to the NL, where enhanced symmetry helps us to 
predict various exact results, is more likely to lead to reliable results. 

4-2. Duality on the NL 

Consider a specific choice of pi, 

e v Kp (i) 

Vi{K P ) = ! (fc) , (67) 
Z^fc=o e 

which enables us to obtain exact results using gauge symmetry under the NL condition 

K = K p . (68) 

The above set {pi} on the NL is denoted by {p NL }- As in the previous section, we try to 
identify the relation between two transition points on the NL using the condition (J61|) . 
Let therefore Xo NL (-^0> Xo NL (-^0 denote the variables on the NL, 

9-1 

X ^(K) = X o({p^},K) = Y,Pi{K)e nV ^ l \ 

1=0 

q-l /<Z-1 \ n 

Xo NL W = XU{P^},K) = q^ 2 J2Pi(K) [J2e v « (k+l) . (69) 

1=0 \k=0 J 

Our conjecture for the replicated systems is the following. If uniqueness of the critical 
point is assumed, the location of it is determined by 

Xo NL (K) = X*o NL (K). (70) 
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If there are two critical points Ki, K 2 on the NL, they satisfy 

Xo NL (^i)Xo NL (^) = Xo NL (^i)Xo NL (^ 2 ). (71) 

We expect these relations to hold even in the quenched (n — > 0) limit. In this limit, 
the condition reads, corresponding to equation (JZDJ), 

g-i , 
-J2pi(K)\og gPl (K) = -, (72) 

1=0 

and to equation (J7TJ) 

9-1 

-^( W (^i)log 9 w(^i)+w(^ 2 )log gW (^ 2 )) = 1- (73) 

It is interesting to note that these conditions are written in terms of the entropy function. 
This fact is suggestive of some underlying structure behind the scene, but we do not 
have a clear interpretation of this fact at the moment. 
It is convenient to define 

oo 

F K (x) = e-T^- 2 ™) 3 (74) 

to facilitate more compact expressions for duality relation. Clearly, we have 

e vx« = jr* (— V (75) 



Q 



We also see 



k=0 k=0 meZ k=0 ^^ 7r ^ y 



q 



e~^ = Yl e ~^~ = F */AQ), (76) 



so that we have 



= F X ly (77) 

Hence equations (J7nj) - (J73|) can be rewritten in terms of Fk(x). It should be noticed 
that the function Fk(x) can be written in terms of $3, a Theta function: 

F K (x) = e~^ x2 Y^ie-^^r 2 ^)^ = e-^ x2 ^Q m2 z 2m = e~^-& 3 (u,Q), (78) 

mgZ m£Z 

where z = e nKx = e m (u = —inKx) and Q = e~ 2n2R . This expression is useful for 
numerical evaluations of the transition points on the computer with special functions 
software implemented. 

Now we would like to discuss the plausibility of our conjectures ()70|) - (J73J). The 
q = 2 and the q = 3 cases are equivalent to the Ising model and the three state 
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Potts model, respectively, both of which have been discussed in detail in These are 
favourable facts to support our conjectures. Unfortunate aspect about the Z q model is 
the difficultly of the analysis of the replicated model even for n = 1. Hence we do not 
give such evidence as we did for the Z 2 case. Instead, the limit q ^> 1 is discussed in 
the next subsection. 



4-3. The limit q — >• oo 

When q is large and the transition point is unique, it is reasonably expected that the 
transition point is of order 0(g). In fact, if we set K = jq and suppose q is large, we find 
a consistent solution as follows. The functions appearing in (J77j) behave asymptotically 

as 

wo) - E> * - ■ <™> 

F 7 , = J2 e-i^- 2 '^* ~ e-f»> 2 . (80) 

Substitution of these expressions into (J72|l leads to 

^^S e ~ uw ^ m - (8i) 

If we approximate the summation by an integral, which should be valid for large q, we 
find 

/OO 1 
dyj(2ny) 2 e-^ 2 = -, (82) 

where we set y — 1/ Jq. We therefore have 

log(2vr 7 ) = 1, (83) 

implying 7 = e/2ir with e being the base of natural logarithm. Hence, when uniqueness 
of the critical point is assumed, the asymptotic location is 

K = — . (84) 
2tt v ; 

When there are two critical points K±,K2, the above argument should be modified 

to some extent. If we assume that both of K\ and K2 are of order O(q), we may set 

Ki = jiq (i = 1, 2) and apply the same procedure as above. The result is that K\ and 

K2 should be related through 

K,K 2 = (85) 

However, for the non-random case, it is known that the two transition points are not 
of order 0(g); one is of order 0(q 2 ) and the other 0(1) [T7] . Since this is expected to 
persist for the random case (see subsection l4.4l below). equation may not necessarily 
capture the true asymptotics of the transition points. 
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4-4- Structure of the phase diagram 

In this subsection we discuss the phase structure of the random Z q model. For the case of 
the non-random model, it was shown that, when q is sufficiently large, the simple Ising- 
like two phase picture (with ferromagnetic and paramagnetic phases) is not possible 
[T7j . The authors of [Hj showed that the transition point determined by the duality 
in the case of only two phases is inconsistent with a kind of Griffith's inequality. The 
discussions were based on the assumption of the existence of the KT transition in the 
continuous model q — > oo. 

We would like to address the issue of the full phase structure of the random Z q model 
but let us restrict our main interest on the NL for a moment. Basically it is expected 
that the phase transitions on the NL are of a similar nature to the non-random case. 
There are three phases; disordered, KT and ferromagnetic. This conclusion may be 
drawn by using the same arguments as for the non-random case. There are, however, 
still some debates about the existence of the KT phase when randomness is introduced in 
the model [23 I2D E2] ■ Hence, in this section, we give a different argument supporting 
the three-phase picture without assuming the existence of the KT transition in the 
continuous model. A crucial point in our arguments is the proof of existence and non- 
existence of the ferromagnetic phase in appropriate parameter regions when q ^> 1. 

First we show the existence of a ferromagnetic phase near the limit K, K p — ► oo 
(the ground state of the non-random system) following [23 [21] • Let us first consider 
the non-random case. The order parameter is bounded as 



(e i 



~)k = Y^P k e 2 < =Po+J2 P ^ m ~ q * P ° " E Pk = 1 - 2 ^ " P °) 

fc=0 fe(^O) fc(^O) 

= 1 - 2(N') K /N, (86) 

where ( )k means thermal average, Pk is the probability that the spin takes the value 
k on each site, and N' is the number of sites such that k ^ 0. All boundary spins have 
k = 0. As in |2Hl|2lj, N' is bounded as 

E (!) 5>. w . < 87 > 

6=4,6,- ^ ' j=l 

where xff is 1 if the jth border (separating a domain of sites from those with different 
fc's) of length b occurs in a configuration and otherwise. The symbol v(b) stands for 
the number of possible borders of length b. The average (X^)k is bounded by, for large 
<?, 

(X^) K < e--(v) 2 . (88) 

Using 

v(b) < 4-3 b q b N/(3b), (89) 
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where the factor q b comes from the possible number of boundary bonds, we find 

< 1 E b3Ve 



N 



Kb I 2tv ' 
2 V q , 



(90) 



6=4,1 



Hence if we choose k = 3ge 2 v <? / sufficiently small, the order parameter, the left-hand 
side of fJ8f)jl . is certainly positive. From this follows that the ferromagnetic phase exists 
when T < C[ 0) /q 2 with C{ 0) being some positive constant. 

Following j23], we can generalize the above argument to the random case. At least 
for K and K p very large, it is possible to prove the existence of a ferromagnetic phase. 
Let us take a border of spin configurations such that X^' ^ and notice that the 
thermal average (X^)k is written as 



Sin n 



{k}LL(ij) 



?V K {ki-kj-lij) 



(91) 



E{fc} ri(y> 

where the sum in the numerator is taken over states in which the jth border of 
length b occurs. Here Uj is the quenched randomness and occurs with probability r 1 ^ 
(r = e~ 2n2Kp / q2 ) approximately at each bond. Take b = 2n (n to be distinguished from 
the replica number) and consider the case where the number of such bonds with l^j ^ 
is m. When m = n,n + l,---,2n, one restricts the summation in the denominator 
of ()91|) to configurations which appear in the numerator and obtains a trivial upper 
bound {X^n) < 1- On the other hand, when m — 0, 1, • • • , n — 1, one finds an upper 
bound pQn) — A 2ri_m_ ^^ with A = e - 27T2K /i 2 ^ where the summation in the exponent 
is over the bonds with l^j ^ 0. This is obtained by restricting the summation in the 
denominator in (jfJlj) to configurations with all ki = 0. Hence the configurational average 



[(^2n)if]av is bounded as 



n-1 
m=0 



2n 



m 



2n 

J2r^ l hX 2n ~ m -^ l h + J2 

lij m=n 

When K < K p , one has r/A < 1 so that this can be replaced by 

2n 



2n 



m 



q m r m . 



(92) 



[<*£?>*u < E 



m=0 



2n 



q m f m \ 2n ~ 2m 



2n 



??? 



g m r m . 



(93) 



Then, if we choose K p /q 2 > K/q 2 ^> 1, by almost the same argument as in the non- 
random case, (N')k/N can be made sufficiently small so that the order parameter does 
not vanish. In particular, there exists a positive constant C\ such that the order is 
nonzero when T < C\jq 2 on the NL. The ferromagnetic phase is also expected when 
Kp/q 2 > K/q 2 ^> 1 but we need a little tighter estimation to prove it. 

On the other hand, to show non-existence of ferromagnetic phase at high 
temperature, we use the arguments in J2H]. We show below 

















(e " )k 








av 



< (e 



2nikj/q\ j 2nikj/q\ 

I if p ,nonrandom \" / i-S^nonrandom; 



(94) 
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where the average on the right-hand side is taken for the non-random model at inverse 
temperatures K p and K, respectively. To verify this inequality, we first notice 



k ■ 




1 2ni-2-\ 




(e o ) K 






av 



gVff^fejj kj 2 lj 1 .j 2 ) 




P V Kp ihx ,32 



Tr 



W e 



2nikj/q 



n 



0'lj2> 



Applying the gauge transformation 



(n>i2 — (71J2 — K ji + K . 



Tr W n - 1J2> e 



VK\kj 1 kj 2 lj 1 j 2 ) 



n i 



we find 



/ 2iri^-\ 

\ e q )k 



(e 9 



k ■ 

)K p {e «)k, 



(95) 

(96) 
(97) 

(98) 



where denotes a certain probability distribution. If we apply a very reasonable 
inequality that the value of the order parameter of the random system is smaller than 
its non-random counterpart with the same inverse temperature 



(e 2wik J /q ) K \ 



^ (c ^ ) K" p ,nonrandom ifi )_R",nonrandom (99) 

to each term in ()98j1 . we arrive at ()94j) . Unfortunately we have not succeeded in proving 
the inequality (|9*9^) mathematically except for the g = 2 case, which was already proved 
in [2E]. However, even in the absence of a formal mathematical proof, equation 
should certainly be valid. Now, if we let K\ denote the phase transition point of the 
non-random model at which the order vanishes, we see that the ferromagnetic phase 
cannot exist when K < K\ or K p < K\. In fact K\ is known to be of order 0(q 2 ) 
and we conclude that the ferromagnetic phase cannot exist on the NL when T > CVg 2 
with C*2 some positive constant. 

Notice that the above arguments imply the existence of (at least one) transition 
point (s) between ferromagnetic and non- ferromagnetic phases. Now if there is a unique 
transition point, it is of order 0(l/q) [cf (|84j)]. But this contradicts with the fact that 
the order vanishes when T > C^/g 2 . Therefore we have shown that the simple Ising- 
like two phase picture is not possible when q is large enough. This conclusion strongly 
suggests that there are three phases in the full phase diagram. A schematic phase 
diagram expected from the above arguments is given in figure 3. 

As a remark, for q not very large, the situation may be more subtle. Numerical 
simulation results suggest that something peculiar may be happening for these models 
(e.g. NER analysis in ^7\). Further careful analyses are necessary on this point. 



2-Kikj/q\ 



4-5. Maximally Random Case 

Let us consider the maximally random case (K p = 0), in which the possibility of a finite 
temperature phase transition has been discussed in the literature [2E|- Here we give an 
argument against such a possibility using duality. 
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oo 



Figure 3. Phase diagram of the random Z q model. Solid lines show expected 
phase boundaries. 



As in the case of the Ising model |3j , the n-replicated system with K p = is related 
to the (n — l)-replicated system on the NL. More precisely, it is easy to prove 

Z n (K, K p = 0)(x Z n ^(K, K P = K), (100) 

where Z n (K, K p ) stands for the averaged partition function (jHJ). Now let us assume the 
applicability of this relation down to the n — > limit and discuss the phase transitions 
for the quenched model with K p = 0. Note that the n — > limit for K p = is equivalent 
to taking the n — > — 1 limit on the NL according to (jlOOj) . 

If we assume that there is a unique transition point, it is determined by the relation 
(jTHj) . which reads 

sr^q-l p nV K (l) 

= q- n+1 . (101) 

Let us take the n —>■ limit. If we set K = c/n, we see 

'9-1 \ n 

J2 eVK(l) ) ~ eVc{0) > ( 102 ) 

J=0 / 

so that (jlOlj) becomes 

This is nothing but a condition to determine the transition point for the non-random 
system when there is a unique transition point. Accordingly we know that c = q/(2n) 
by (|66j) with K\ = K 2 and hence K = c/n — > oo as n — ► 0. This means that the phase 
transition, if it is unique, occurs at T = 0. The same discussions can be applied to the 
case where there are two transition points. In any case the conclusion is that the phase 
transitions, if there exist, occur only at T = 0. 
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We should remember that our discussions were based on the subtle assumption 
about the applicability of the relation as n — > (n — ► — 1 on the NL) and the assumption 
that (J70|) or (J7T|) gives the criticality condition. Hence our conclusion has not been 
completely rigorously derived. It would be interesting to clarify this point with other 
methods or arguments. 

5. Conclusion and discussion 

In this paper, we have generalized the duality argument combined with the replica 
method, which was originally applied to the 2c? random bond Ising model on the square 
lattice, to a variety of random spin systems. Our main results are the conjectures on 
the transition points of systems with quenched randomness. 

First we considered the random Z2 models. We have given the conjecture (|21|) 
about the relation between the multicritical points of two models which are mutually 
dual. Besides exact computations for the replicated system with n = 1, 2 and 00, 
numerical simulations support our conjecture for the 2d ±J random bond Ising model 
on the triangular/hexagonal lattices, the 3d random Ising/gauge models on the cubic 
lattice and the 2d anisotropic random bond Ising model on the square lattice. We think 
it remarkable that a single theoretical framework makes it possible to derive a series of 
predictions to be compared favourably with many independent numerical simulations. 

Next we treated the random Z q model. We have shown that there are at least two 
transition points for sufficiently large q. The most probable scenario is that there are 
three phases in the phase diagram, paramagnetic, KT-like and ferromagnetic phases. 
Our arguments, however, do not assume the existence or some specific properties of the 
KT phase. By applying the duality argument to this model, we have also given the 
conjectured relation (|73j) between the two transition points. 

An interesting question is why we restrict ourselves to the NL. The fixed-point 
condition of the leading Boltzmann factor xq(K, p) = Xq(K,p) relates p and K and 
may give the whole shape of the phase boundary. This is indeed the case for n = 1 
as well as for n = 2 above the multicritical point. We nevertheless have restricted 
ourselves to the NL in this paper because it is difficult to directly verify the ansatz 
on the whole part of the phase diagram numerically for many systems with quenched 
randomness. It should also be kept in mind that the multicritical point is the place 
where two completely different types of symmetries, invariance under duality and gauge 
transformation, meet under the present conjecture. This implies that the multicritical 
point has clearly distinguished symmetry features, which allows us to discuss this point 
on a different basis from other points of phase transition. 

We believe that sufficient analytical and numerical evidence has been accumulated 
to support the validity of our conjecture. It is an interesting future problem to provide 
a mathematically rigorous proof. 
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Appendix 

In this appendix we calculate the prefactor in the duality relation of the partition 
function (jSJ) in section 2. Here we follow the derivation of the duality relation by Wu 
and Wang [Bj. 

The partition function of the system described by the Hamiltonian (PJ in section 2 

is 

Z = E W U ^c\ (104) 

where 

u&c) = II ex P &0 • ( 105 ) 

xedc 

U(C,c) is the Boltzmann factor for an element C with dimension r. For example, the 
usual two-body interaction on a bond has r = 1 and the lattice gauge theory has r = 2. 
£x is a modulo-2 spin variable which takes or 1, and £c* is defined by £c* = ^2xeac^ x 
modulo 2. 

Let us define the dual Boltzmann factor for the dual element C* by the binary 
Fourier transformation, 

U*(\ c *) = 2^ exp(27r^ c Ac7*M£c), (106) 
Cc=o,i 

or conversely, 

U{ic)=2^ exp(27riA c *£ c )[/*(A c *) 

A C *=0,1 

= 2-3 exppTriAo ^ Zx)u*(*c*) 

a c *=o,i V xeac J 

= E ( II T(&,\c*))u*(\c*), (107) 
a c *=o,i \xedc J 

where 

T{£ x , Ac*) = T^k exp(2vrzAc^ a; ), (108) 
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and Be is the number of (r — l)-dimensional elements on the boundary of C. 

Inserting (jlU7|) into (jlU4|) and taking the sum over we obtain a modulo-2 
Kronecker delta for each 

£ n ntx,*c) 

= £ II 2"^exp(2«A c *^) 
£25=0,1 c* : a;e9c 

= 2 i-Ec:te C 4 Ld2 | ^ Ac ,j (109) 

It is useful to define the following symbol of a constrained sum, which stems from the 
Kronecker deltas in 



£ = £ n^M £ A H- ( n °) 

{A C *=0,1} {A C »=0,1} X \C*:X£dC 



Using this, (jl(J4|) becomes 



z= ( ii^^^k) e' ) ■ ( m ) 

a; / V {a c ,=o,i} c 



The prefactor can be simplified as 

2 1 ~Sc : £Eeec 

x 

_ 2 Nr ~ 1 ~k^'X T-c-.Xeac ^ 

= 2*- 1 -^, (112) 



where the number of m-dimensional elements is denoted by N m . 
The final result is 

z = 2 N r -,-^ Y[U*(\ C *). (113) 

{A C *=0,1}C* 

Here we have replaced the product over C with C*, which are identical operations. 

The partition function (jlU4|) can also be rewritten using the Boltzmann factors for 
the element C, 

Z = E U U ^) (H4) 
{€a;=o,i} c 

=2^ e' n^)- ( iis ) 

{?c=0,l} c 

Note that U(£c) i n (|H4|) is a function of £x as in (1), but it is not in (|115j) : It is 
redefined on each element C as a function of £c itself. The usual sum over £ in ()114|) 
can be written as in f)l 15|) using the definition (jllOj) if we consider a correct mapping 
from the configuration space defined by onto the one by £c- The prefactor 2^ 9 in 
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(|115|) is the degree of the ground-state degeneracy. The mapping from the £x space to 
the £c space is not one-to-one. In (jll5J) the ground state configuration is £c — for 
all C for ferromagnetic interactions and is unique, while in (jll4j) the ground state is 
degenerate. Thus the mapping is 2 A/ ' 9 to I. For example, 2 A/ ' 9 = 2 for r = 1 because 
all-up and all-down states are degenerate. For r = 2, 2^0 is dependent on the number 
of sites because the model has Z 2 gauge symmetry in this case. Wegner calculated this 
degree of degeneracy in [5] under general conditions and the result is 

r-2 

Af g = ^(-l) r " m iV m + ti for r > 2, 

m=0 

= 1 for r = 1, (116) 

where t\ is a constant which depends on the topology of the lattice (or the boundary 
condition) and not on the number of elements. If we use the generalized Euler's relation 
for the number of lattice elements [5], 

d 

J2(-l) m Nm = t 2 , (117) 

m=0 

where t% is also a constant dependent on the topology, M g becomes 

d 

K= E (-^^Nm + t, (118) 
m=r— 1 

with t being another constant determined by t\ : t 2 and r. 

We can derive the factor ()118|) intuitively from the difference of the numbers of 
spin configurations between the expressions ()114j) and ([115|) . In (I114|) the number of 
independent £a;'s is iV r _i. In flll5|l the number of independent £c's may appear to be 
N r , but we must take the number of constraints in ()115|) into account, which is iV r+1 
and should be subtracted from iV r . However these constraints are redundant or not 
independent of each other because r + 1 dimensional element is always on the boundary 
of r + 2 dimensional element. Therefore we must subtract N r+ 2 from the number of 
constraints AT r+1 . However the number N r+ 2 is also redundant again and N r+ 3 must 
be subtracted from it, which is similar for higher dimensional elements. Then M g is 
calculated as 

M g = N r ^ - {N r - (N r+1 - (N r+2 - (N r+3 . . .} 
d 

= hl) m - {r - 1} N m + t, (119) 

m=i — 1 

which is equivalent to (|118|) . The constant t will be the same as in (|118|) . 

Using (jllHj) . pi5j) and ()118jl . we obtain the duality relation between the original 
and the dual partition functions, 

2^ Ylll U &) (=Zor is {MK),u^(K)}) 
Uc=o,i} c 

= 2 N g * +a ]JU*(\c*) {=2 a Z dnal {ul(K),u*_ l (K)}), (120) 

{A C *=0,1}C* 



a = -M; + N r _i - -jf . (122) 
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where 

r+l 

^ = ^(-l) m_(r+1) iV m + t*, (121) 

m=0 

and t* is also a constant which depends on the topology, a is defined by 

2 

Remember that we use the modulo-2 spin variables here, while spins take the value ±1 
in the main text. 

Now let us consider the self-dual case. In this case d must be even and r = d/2. 
Furthermore the number of m dimensional elements satisfies 

N m = N d _ m . (123) 

Inserting this into the Euler's relation ()117j) . we obtain 

d/2_1 AT d N + 

J2(- 1 ) m N m + (-l) d/2 ^= £ (-l) m N m +(-l) d/2 ^l = % (124) 

m=0 m=d/2+l 

In addition, Af g and Af* become 

M 9 = £ (-^ W2 ^N m + t = N d/2 _, + (-lp- 1 ! + t, 

m=d/2-l 



d/2+1 



a/"* = £ (-i) m -w 2+1 )iv m + r = Ar d/2+1 - + (-i) d / 2+1 | + f 



Hence, 



= N d/2 ^ - ^ + (-I)' 72 - 1 ! + f = A^ 9 + 1* - 1. (125) 

a= (-l) d/2 | -t*. (126) 

Therefore a is just a trivial constant in the self-dual case and negligible in the 
thermodynamic limit. The factors 2^ 9 and 2 A/ 9 in ()120|) . which differ only by a 
trivial constant, do not concern when we derive the transition point from the relation 
u±i(K c )=u* ±1 (K c ). 
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